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ABSTRACT 

We show that the irradiated accretion disk model can account for all the op- 
tical and infrared observations of the anomalous X-ray pulsars in the persistent 
state. Model fits do not constrain the outer disk radii, while placing an upper 
limit to the inner disk radii, and thus to the strength of the dipole component of 
the stellar magnetic field. While magnetar fields (B* > 10 14 G) in higher multi- 
poles are compatible with the irradiated disk model, magnetic dipole components 
of magnetar strength are not consistent with optical data. 

Subject headings: pulsars: individual (AXPs) — stars: neutron — X-rays: bursts 
— accretion, accretion disks 

1. Introduction 

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray Repeaters (SGRs) constitute 
a special class of neutron star systems (Mereghetti et al. 2002; Hurley 2000; Woods & 
Thompson 2004). They are identified mainly through their X-ray luminosities (L x ~ 10 34 — 
10 36 erg s^ 1 ), which are orders of magnitude higher than their rotational powers Erot — 
Their spin periods are clustered to a very narrow range (5 — 12 s). All five known SGRs 
and two of the eight known AXPs show repetitive, short (^ls) super-Eddington bursts with 
luminosities up to 10 42 erg s _1 . Three giant flares with peak luminosities L p > 10 44 erg s -1 
and durations of a few minutes were observed from three different SGRs (Mazets et al. 1979; 
Mazets et al. 1999; Hurley et al. 1999; Palmer et al. 2005). 

The short time scales and the super-Eddington luminosities of the soft gamma-ray bursts 
strongly indicate a magnetar mechanism for these bursts. Magnetar models (Duncan & 
Thompson 1992; Thompson & Duncan 1995; Thompson & Duncan 1996) adopt strong 
magnetic fields with magnitude > 10 14 G on the stellar surface to explain the burst 
energetics. Are the burst energies stored in the dipole component or the higher multipoles of 
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the magnetic field of the neutron star? In the current magnetar models, the dipole component 
of the magnetic field must be of magnetar strength to account for the spin down properties 
of the AXPs and SGRs. In these models, persistent X-ray luminosities are explained by the 
magnetic field decay, while the magnetic dipole torque is taken as the mechanism responsible 
for the spin down rates of the sources. 

In the alternative fallback disk model (Chatterjee, Hernquist, & Narayan 2000; Alpar 
2001), the source of the X-rays is accretion onto the neutron star, while the optical/IR light 
originates from the accretion disk. Rotational evolution of the neutron star is determined 
by the interaction between the disk and the magnetosphere of the neutron star (B* ~ 
10 12 — 10 13 G). Fallback disk models can account for the period clustering of AXPs and 
SGRs as the natural outcome of disk-magnetosphere interaction during their lifetimes (Alpar 
2001; Ek§i & Alpar 2003). These models are consistent with magnetar fields on the neutron 
star provided that these fields are in higher multipole components of the magnetic field. As 
higher multipole fields rapidly decrease with increasing radial distance (as r -5 for quadrupole 
component), it is the dipole component of the magnetic field which determines the interaction 
and the angular momentum transfer between the disk and the neutron star. The strength of 
magnetic dipole field in order to explain the period clustering of the AXPs and SGRs over 
their M history is 5* ~ 10 12 - 10 13 G (Alpar 2001; Eksi & Alpar 2003). 

Fallback disk models can also explain the enhancements observed in the persistent lu- 
minosities of SGRs and AXPs. The X-ray enhancement of the SGR 1900+14 following its 
giant flare can be explained by the relaxation of a disk which has been pushed back by a 
preceding burst (Ertan & Alpar 2003). The same model with similar disk parameters can 
also reproduce the correlated X-ray and IR enhancement of AXP 2259+58, which lasted for 
~ 1.5 years, if this is triggered by a burst, with a burst energy estimated to have remained 
under the detection limits (Ertan, G6gu§ & Alpar 2006). 

The suggestion of fallback disks has motivated observational searches for disk emission 
in the optical and IR bands, and resulted in various constraints on the models. Some of the 
AXPs were observed in more than one IR band (Hulleman et al. 2001; Israel et al. 2002; 
Wang & Chakrabarty 2002; Kaspi et al. 2003; Israel et al. 2003; Hulleman et al. 2004; 
Israel et al. 2004; Tam et al. 2004; Morii et al. 2005; Durant & van Kerkwijk 2006b). 
AXP 4U 0142+61 is the source with most extended observations, as it was also observed 
in the optical R and V bands (Hulleman et al. 2000; Hulleman et al. 2004; Dhillon et al. 
2005), and recently in mid-IR bands with the SPITZER observatory (Wang et al. 2006). 
Discovery of modulation in the R band luminosity of 4U 0142+61 at the neutron star's 
rotation period P=8.7 s, with a pulsed fraction 27 % (Kern & Mertin 2002; Dhillon et al. 
2005), is particularly significant. This fraction is much higher than the pulsed fraction of the 
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X-ray luminosity of this source, indicating that the origin of the pulsed optical emission is not 
likely to be the reprocessed X-rays by the disk. Magnetospheric models for these pulsations 
can be built with either a dipole magnetar field or within a disk-star dynamo model (Cheng 
& Ruderman 1991), in which a magnetospheric pulsar activity is sustained by a stellar dipole 
field of ~ 10 12 G and a disk protruding within the magnetosphere. Ertan & Cheng (2004) 
showed that this pulsed optical component of the AXP 4U 0142+61 can be explained by both 
types of magnetospheric models. Thus the presence of strong optical pulsations from the 
magnetosphere does not rule out the possibility of a fallback disk together with 10 12 — 10 13 
G surface dipole magnetic field. 

In the present work, we concentrate on the unpulsed optical/IR emission from the AXPs 
and SGRs in their persistent states, and test the expectations of the irradiated accretion 
disk model through the observations in different optical/IR energy bands (V, R, I, J, H, 
K, and K s ). The optical/IR emission expected from the irradiated fallback disks was first 
computed and discussed by Perna, Hernquist & Narayan (2000) and Hulleman, van Kerkwijk 
& Kulkarni (2000). Using similar irradiation strengths, Perna et al. (2000) and Hulleman 
et al. (2000) found similar optical fluxes that remain well beyond those indicated by the 
observations of AXP 4U 0142+61 and AXP IE 2259+586. To explain this result, Perna et 
al. (2000) suggested that the inner disk regions could be cut by an advection dominated 
flow, while Hulleman et al. (2000) concluded that the then existing optical data of the AXP 
4U 0142+61 (in I,R,V bands) can only be accounted for by an extremely small outer disk 
radius, around a few xlO 9 cm. In the present work, we show that the optical/IR data of 
the AXPs can be explained by the irradiated accretion disk model without any implausible 
constraints on the outer and inner disk radii. The main reason for the difference between our 
results and those of earlier works is that both Hulleman et al. (2000) and Perna et al. (2000) 
assumed a particular irradiation strength, while we keep it as a free parameter, to address 
the broad band full data set. This approach is supported by the observations of the low 
mass X-ray binaries (LMXBs) which indicate varying irradiation strengths. Furthermore, 
model fits are sensitive to the interstellar reddening parameter Ay, which was estimated to 
be between 2.6 and 5.1 for 4U 0142+61 (Hulleman et al. 2004). For this source, we obtain 
the best model fit with Ay = 3.5 which turned out to be consistent with the recent result 
Ay = 3.5 ± 0.4 (Durant & van Kerkwijk 2006a). On the other hand, the test disk model 
with unreasonably small outer disk radius requires Ay = 5.4 (Hulleman et al. 2000). In this 
model, optical emission comes from the outer disk, while in our model it is the inner regions 
of an extended disk that emits substantially through the optical bands (see § 3 for further 
discussion). We give the details of the disk model in § 2. We discuss our results in § 3, and 
summarize the conclusions in § 4. 
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2. Optical/IR Emission from the Irradiated Disk 

Model fits to the X-ray and IR enhancement data (Kaspi et al. 2003) of AXP IE 
2259+586 favor the irradiated disk model, though they do not exclude the nonirradiated 
thin disk model (Ertan, G6gii§ & Alpar 2006). We start by assuming that the AXP disks are 
irradiated and include the irradiation strength as a free parameter through our calculations. 

When the disk is irradiated by the X-rays from the neutron star, both the intrinsic 
dissipation and the irradiation flux should be taken into account in calculations of the disk 
blackbody emission. A steady disk model is a good approximation for the present evolution 
of the AXP and SGR disks in their persistent states. For a steady thin disk, the intrinsic 
dissipation can be written as D = (3/8tt)(GMM/R 3 ) (see e.g. Frank et al. 2002) where M 
is the disk mass flow rate, M is the mass of the neutron star and R is the radial distance 
from the neutron star. In the absence of irradiation, the effective temperature T e g of the 
disk is proportional to R~ 3 ^ for a given M. For an irradiated disk, the irradiation flux can 
be written as F- irr = uT^ T = (CMc 2 )/(AttR 2 ). (Shakura & Sunyaev 1973), where c is the 
speed of light. Iraadiation parameter C includes the effects of the conversion efficiency of the 
accretion into X-rays, disk geometry and the albedo of the disk face. Irradiation temperature 
T irr = (Firr/cr) 1 / 4 is proportional to R~ 1 / 2 . For small radii, dissipation is the dominant 
source of the disk emission. At a critical radius R c , the irradiation flux becomes equal to the 
dissipation rate, and beyond R c , the disk emission is supported mainly by reprocessed X-rays. 
Equating F irr to D, the critical radius is found to be R c = (3GM,)/(2Cc 2 ) ~ (10- 4 /C)3 x 10 9 
cm. The effective temperature profile of the disk can be obtained using aT 4 ff = D + -Firr 
where a is the Stefan-Boltzmann constant. 

We adopt the observed magnitudes in the optical/IR bands, distances and the N H values 
given by Woods & Thompson (2004) and references therein and convert the magnitudes to 
energy flux values. We calculate Ay values using Nh = 1-79 x 10 21 Ay (Predehl & Schmitt 
1995). To find the model disk flux in a given observational band, we integrate the calculated 
blackbody emissions of all radial grids radiating in this band. For comparison with data, we 
calculate the model disk fluxes along the optical/IR bands V, R, H, I, J, K and K s . For all 
sources we set cos % — 1 where i is the angle between the disk normal and the line of sight of 
the observer. We equate the disk mass flow rate M to the accretion rate onto the neutron 
star, thus assuming the mass loss due to the propeller effect is negligible. We first adjust 
M to obtain the observed X-ray flux. Next, using this value of M and taking the strength 
of the magnetic dipole field = 10 12 G on the surface of the neutron star we calculate the 
Alfven radius Ra which we take to be the inner radius of the disk. Then, we look for a good 
fit to the overall available optical/IR data by adjusting the irradiation strength C within the 
uncertainties discussed in 8 3. 
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3. Results and Discussion 

Our results are summarized in Table 1. For each source, the first column gives the 
unabsorbed flux data obtained from the observed magnitudes and the estimated Ay values 
(see Table 1) given in Woods & Thompson (2004), and the second column gives the model 
fluxes. For the AXP 4U 0142+61, the range of reddening quoted in earlier literature is 
2.6 < Ay < 5.1 (Hulleman et al. 2004). We obtain a good fit with Ay = 3.5. Table 1 shows 
that the irradiated steady disk model is in agreement with all the AXPs observed in the 
optical and IR bands. The parameters of the model for each source are given in Table 2. 

At present, AXP 4U 0142+61, which has been observed in five different optical/IR 
bands from K to V in the same X-ray luminosity regime, seems to be the best source to 
study the properties of AXPs in the persistent state. Earlier work by Hulleman et al. (2000) 
excluded the disk model for the AXP 4U 0142+61. They obtained an irradiation temperature 
profile by using a particular irradiation strength. The estimated optical flux for an extended 
disk with this irradiation efficiency remains above the optical data points of the AXP 4U 
0142+61 (see Fig. 3 in Hulleman et al. 2000). Considering the possibility that the optical 
flux might originate from the outermost disk region, Hulleman et al. (2000) tried to fit the 
then observed three data points in the I, R and V bands to the Ray leigh- Jeans tail of a 
blackbody spectrum with the extinction parameter Ay = 5.4. This placed an upper limit to 
the outer disk radius which is too small for a realistic disk. The key factor in the difference 
between the earlier results and our recent results is the irradiation efficiency, which we allow 
to vary in conjunction with Ay, to provide the best fit to the current broad band data. 
We note that the irradiation efficiency indicated by the observations of the low mass X-ray 
binaries varies from source to source. Even for the same source, the ratio of the irradiation 
flux to the X-ray flux may change with accretion rate (de Jong et al. 1996; Dubus et al. 
1999; Ertan & Alpar 2002). Taking these into account, we keep the irradiation efficiency as 
a free parameter for our model fits. With the parameters given in Table 2, the irradiated 
disk model can account for the optical/IR data of this source without setting any stringent 
constraints on the inner or outer disk radii. In our model, the optical luminosity is radiated 
from the inner disk, while longer wavelength IR emission comes from larger radii. A more 
detailed analysis of the AXP 4U 0142+61 with the new detections in the mid-IR SPITZER 
bands confirms the results here (Ertan et al. 2006). The irradiation parameter C obtained 
from our model fits turned out to be in the range (10 -4 < C < 10~ 3 ) estimated from the 
observations of LMXBs and the disk stability analyses of the soft X-ray transients (de Jong 
et al. 1996; Dubus et al. 1999; Ertan & Alpar 2002). Within the critical radius R c given by 
Eq. 4, dissipation is the dominant heating mechanism. For the disk model of the AXP 4U 
0142+61, R c ~ 3 x 10 cm and R m = 1 x 10 9 cm. The innermost disk emitting mostly in 
the UV bands also contributes to the optical emission. The radial distance at which the disk 



- 6- 



blackbody temperatures peak at the optical bands (R,V) is about 10 10 cm. at this radial 
distance, about 35 % of of the optical radiation is due to dissipation, and the rets is due to 
irradiation. Peak temperatures of the IR bands from I to K s lie between R ~ 2 x 10 10 cm 
and R ~ 1.5 x 10 11 cm. 

There are several uncertainties related to the inner disk emission characteristics of the 
AXPs, which are not possible to address by the irradiated thin disk model. Fistly, emission 
properties of the innermost disk boundary interacting with the magnetosphere are not very 
clear. Secondly, the contributions from the magnetospheric pulsed emission which is known 
to have a fraction about 27 % in the R band for 4U 0142+61, is likely to be radiated from 
the other IR and optical bands as well. Relative amplitudes of these pulsed contributions 
radiated from different optical/IR bands are not known at present. Finally, there could be 
some X-ray shielding effects depending on the details of the geometry of the innermost disk 
regions, which could also affect the optical/IR emission properties of these sources. For all 
the AXPs that were detected in the optical/IR bands, optical and IR flux values of our 
models remain within about 30 % of all the data points, which is a reasonable fit considering 
the uncertainties discussed above. 

For AXP J 1708-40, Durant and van Kerkwijk (2006b) recently found that the previously 
reported IR data in K s , H, and J bands are likely to be a background star. They found 
another object within the positional error cycle and argued that this second object is more 
likely to be the IR counterpart to the AXP J1708-40. For this source, we adopt the IR (K s , 
H, J) data set reported by Durant and van Kerkwijk (2006b). 

For the persistent state of the AXP IE 2259+586, we use the preenhancement data 
(Hulleman et al. 2001). This source was detected in K s band and there are upper limits 
for I and R bands. Our model flux values are three and ten times below the upper limits 
reported for I and R bands respectively. 

AXP IE 1048-59 was detected in K s , H and I bands (Wang & Chakrabarty 2002). Ob- 
served X-ray flux from this source between December 2000 to January 2003 show a variation 
within a factor of 5 (Mereghetti et al. 2004). We use the X-ray flux obtained from the 
nearest X-ray observation to the date of the IR observations. 

AXP IE 1841 was detected only in the K s band, and there is a high upper limit in the 
R band (Watcher et al. 2004). Model estimates in other optical/IR bands for this source 
(and the other AXPs) can be tested by future optical and IR observations. 

Since there are no detections in short wavelength optical bands for the AXPs (except 
for AXP 4U 0142+61), model fits are not sensitive to the chosen inner disk radii. We equate 
the inner disk radii to the Alfven radii (Table 2) corresponding to a magnetic field with 
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magnitude = 10 12 G on the stellar surface and the accretion rates derived from the 
estimated X-ray luminosities (see Table 1 for references). For the AXP 4U 0142+61, optical 
data in R and V bands provide a constraint for the inner disk radius, and thereby for the 
strength of the magnetic dipole field of this source (see § 4). 

4. Conclusion 

We have shown that the optical, infrared and X-ray observations of the AXPs in their 
persistent states can be explained with irradiated disk models. Among the AXPs, 4U 
0142+61 is currently the only source which provides an upper limit for the inner disk radius 
through its optical (R,V) data. For the best model fit for this source, which we have obtained 
with Ay = 3.5, the model inner disk radius (~ 10 9 cm) is around the Alfven radius for the 
accretion rate, estimated from the X-ray luminosity, together with a dipole magnetic field 
strength B* ~ 10 12 G on the neutron star surface. Nevertheless, it is possible to obtain 
reasonable fits by increasing the inner disk radius and decreasing the reddening accordingly. 
For Ay = 2.6, the minimum value of the reddening in the range 2.6 < Ay < 5.1 (Hulle- 
man et al. 2004), we obtain the best fit with R in ~ 8 x 10 9 cm which corresponds to the 
maximum reasonable dipole field strength B t ~ 4 x 10 13 G on the pole and half of this on 
the equator of the neutron star. We note that these limits could be increased depending on 
the amount of possible mass loss due to propeller effect and/or on how much the inner disk 
radius penetrates inside the Alfven radius (see Ertan et al. (2006) for a detailed discussion 
for 4U 0142+61). On the other hand, even including these possibilities, very recent analysis 
concluding Ay = 3.5 ± 0.4 (Durant & van Kerkwijk 2006a) implies surface dipole magnetic 
field strengths less than about 10 13 G. While the magnetar fields (B* > 10 14 G) in multipoles 
are compatible with this picture, optical (R, V) data excludes a hybrid model involving a 
disk surrounding a magnetar dipole field. In the latter case, inner disk regions emitting in 
the optical would be truncated by the magnetar dipole field. High magnetic fields in multi- 
poles, on the other hand, decrease rapidly with increasing radial distance, and do not affect 
the disk magnetosphere interaction. 

Uncertainties in the source distances, require modification of the model M values, and 
thereby the disk temperature profiles. For the AXP 4U 0142+61, new fits with similar 
quality can be obtained by slightly modifying the model inner disk radius and the irradiation 
parameter C as long as such distance corrections remain within a factor of ~ 3. In this case, 
surface dipole magnetic field should also be recalculated consistently with the newly obtained 
M and Alfven radius, ta oc M' 2 ! 1 B^J 1 . For the other AXPs, which were observed only in 
the IR bands, model fits at present are not sensitive to the inner disk radius, and in the 
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case of distance corrections, similar fits can be obtained only by changing the irradiation 
strength. This is because IR radiation is emitted from the irradiation dominated outer disk 
regions, and therefore M corrections do not change the relative amplitudes of the fluxes in 
different IR bands. Future detections of these sources in the optical bands can constrain the 
inner disk radii together with the surface dipole field strengths. 

On the other hand, existing IR data of the AXPs, including recent observations of 4U 
0142+61 by SPITZER in 4.5 /im and 8 jum bands (Wang et al. 2006), do not put an upper 
limit for the extension of the outer disk radius R ut- The lower limit for R out provided by the 
longest wavelength IR data of the AXP 4U 0142+61 is around 10 12 cm. Further observations 
in the longer wavelength infrared bands by SPITZER space telescope will provide valuable 
information about the structure and possibly the extension of the fallback disks around these 
systems. As a final remark, some AXPs and SGRs which are under the detection limits in 
some of the optical and IR bands could be observed in these bands if they exhibit phases of 
enhanced emission, as observed in the SGR 1900+14 and the AXP 1E2259+586. 
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Note. — The data flux values were calculated by using the magnitudes and Ay values given in the references below. 
For the AXP 4U 0142+61 plausible range for reddening is 2.6< Ay < 5.1 (Hulleman et al. 2004); the data of this source 
here correspond to Ay=3.5. REFERENCES: (J1708-40) Durant & van Kerkwijk 2006b, Rea et al. 2003; (IE 2259+586) 
Hulleman et al. 2001, Woods et al. 2004; (4U 0142+61) Hulleman et al. 2000, 2004, Patel et al. 2003, Morii et al. 2005; 
(IE 1841-45) Wachter et al. 2004, Morii et al. 2003; (IE 1048-59) Wang k Chakrabarty 2002, Mereghetti et al. 2004 
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Table 2: The parameters of the irradiated disk model which gives the optical/IR flux values 
seen in Table 1. For all the sources, we set cos? = 1 where i is the inclination angle between 
the disk normal and the line of sight of the observer, and we take the outer disk radius 
-Rout = 5 x 10 12 cm. See Section 3 for details. 
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